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Abstract 

We perform a microcanonical study of classical lattice ip 4 field models in 3 
dimensions with 0(n ) symmetries. The Hamiltonian flows associated to these 
systems that undergo a second order phase transition in the thermodynamic 
limit are here investigated. The microscopic Hamiltonian dynamics neatly 
reveals the presence of a phase transition through the time averages of con¬ 
ventional thermodynamical observables. Moreover, peculiar behaviors of the 
largest Lyapounov exponents at the transition point are observed. A Rieman- 
nian geometrization of Hamiltonian dynamics is then used to introduce other 
relevant observables, that are measured as functions of both energy density 
and temperature. On the basis of a simple and abstract geometric model, we 
suggest that the apparently singular behaviour of these geometric observables 
might probe a major topological change of the manifolds whose geodesics are 
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the natural motions. 
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I. INTRODUCTION 


The general problem of the relevance of microscopic dynamics to the statistical behaviour 
of physical systems dates back to Boltzmann’s ideas at the very beginning of statistical 
mechanics, and is still far from being clarified and solved. Within this framework, one can 
extract a less general but still challenging question, i.e. whether the microscopic Hamiltonian 
dynamics displays some relevant change when a given system undergoes a phase transition. 

Studying microscopic Hamiltonian dynamics means that - instead of using ensemble sta¬ 
tistical averages - one numerically computes time averages of the relevant observables. There 
are two main reasons for so doing: i) there exist interesting observables that are intrinsi¬ 
cally dynamical, as is the case of Lyapounov exponents; ii) through a differential-geometric 
description of the dynamics, based on simple tools of Riemannian geometry, new concepts 
and methods come to enrich the standard approaches to the study of phase transitions, 
hinting to a possibly deeper characterization of their very nature from the standpoint of the 
mathematical structures involved. 

The geometric formulation of the dynamics of many-degrees-of-freedom systems was first 
used by Krylov in his pioneering studies on the dynamical foundations of statistical mechan¬ 
ics |T[ . Then, during the last two decades, there have been some attempts to cope with the 
ergodicity of Hamiltonian systems through a geometric theory of dynamics P|. A more 
recent series of papers ||-|8|], instead of dealing with ergodicity, successfully addresses the 
problem of explaining and quantifying Hamiltonian chaos within a geometric framework 
where natural motions are seen as geodesics of a suitable Riemannian manifold (henceforth 
referred to as “mechanical manifold”). Here chaotic dynamics stems from curvature fluctu¬ 
ations along the geodesics, through a mechanism similar to the parametric destabilization 
of the stable orbits of a pendulum. At variance with a widespread belief, negative curva¬ 
tures do not appear essential to produce chaos, positive and fluctuating curvatures can work 
as well. A very interesting point is that the average degree of instability of the dynam¬ 
ics is given in terms of curvature-related quantities integrated over the whole “mechanical 
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manifold”. This establishes a link between a dynamical aspect of a given system — the sta¬ 
bility/instability of its trajectories — and some global geometric properties of its associated 
“mechanical manifold”. 

Now, when a model system displays a phase transition, a natural question arises: what 
kind of relationship exists - if any - between all the well-known major thermodynamic 
changes occurring at the transition point and the mentioned global geometric characteristics 
of the “mechanical manifolds”? The present work actually shows that a second order phase 
transition appears to be associated with an abrupt change in the global geometry — and 
possibly in the topology, as we conjecture — of the “mechanical manifolds”. 

The above problematics is addressed in the present work by studying the dynamics of 
classical field theories, discretized on a lattice. A classical lattice field theory can be regarded 
as a discrete classical dynamical system. In particular, we shall consider the classical (p A 
theory, whose lattice version is a set of coupled nonlinear oscillators. 

Equilibrium phase transitions are usually studied in the framework of the Gibbsian 
canonical ensemble. Dynamics, when it is considered, is introduced only a posteriori: the 
most common procedure is to describe it by means of non-deterministic equations, usu¬ 
ally of the Langevin type, whose limiting probability distribution is the Boltzmann weight 
exp(— H/T) where H is the Hamiltonian of the system. Here we are going to adopt a 
completely different approach, i.e. from the very beginning we consider the deterministic 
Hamiltonian dynamics, without making explicit assumptions on the equilibrium properties 
of the system, and we observe how the phase transition is signaled by the dynamics. On 
rigorous grounds, one cannot be sure that a phase transition does exist in a system studied 
through its Hamiltonian dynamics, because there is no proof of the fact that the dynamics 
is ergodic. Moreover, even assuming that it is ergodic, the ergodic measure will be the 
microcanonical rather than the canonical one. The two ensembles are equivalent only in the 
thermodynamic limit, thus the phenomenology observed in finite systems, as the systems 
considered in numerical simulations necessarily are, might be different. To give only an 
example, let us consider the phenomenon of ergodicity breaking, i.e. the fact that ergodicity 
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is not valid for the whole phase space but only for disjoint subsets of it. Such a phenomenon 
is indeed tightly related to phase transitions; in fact, when it occurs it entails a symmetry 
breaking, as in usual phase transitions. But ergodicity breaking is a more general concept 
than symmetry breaking, in fact it is also at the origin of those phase transitions which do 
not correspond to the breaking of an evident symmetry of the Hamiltonian (for example in 
spin glasses) 1*1 [H|. In the canonical ensemble, ergodicity can be broken only in the ther¬ 
modynamic limit |12|, while in the microcanonical ensemble, in principle, there might be 
ergodicity breaking also in finite systems. Ergodicity being a dynamical property, we think 
that a dynamical approach is particularly appropriate to study such a phenomenon. 

It is worth mentioning here that ergodicity breaking in classical Hamiltonian systems 
can be related with supersymmetry breaking 0; this relation is estabilished within the 
framework of a path-integral formulation of classical mechanics, where the bosonic sector of 
a supersymmetric Lagrangian is given by a suitable function of the canonical coordinates, 
obeying standard Hamilton equations, and the fermionic sector contains ghost fields that 
- rather surprisingly - obey the Jacobi equation describing the stability of classical paths 
TJ|.[To|l. In this framework the spontaneous symmetry breaking can occur also at finite 
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Our results show that, as far as the lattice models considered are concerned, the 
numerical phenomenology, obtained by simulating Hamiltonian dynamics, is perfectly con¬ 
sistent with the expectations based on equilibrium statistical mechanics. Moreover, we inves¬ 
tigate whether the instability of dynamical trajectories, measured by Lyapounov exponents, 
is sensitive to the phenomenon of phase transition ||T6[. In the light of the geometrization of 
dynamics, Lyapounov exponents are also seen as probes of the hidden geometry of motion, 
and in fact our results suggest that the deep origin of ergodicity breaking and of the dynam¬ 
ical counterpart of a phase transition could be found in a major change in the geometric - 
or even topologic - structure of the “mechanical manifolds” underlying the dynamics. 

The paper is organized as follows: in Section |IJj we introduce the models studied, we 
describe the numerical techniques which we adopted, and we discuss the phenomenology of 
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the phase transition as it emerges from the dynamics. In Section [ITI] the main definitions and 
results of the Riemannian description of Hamiltonian chaos are given, and the behaviour 
of the geometric observables in our models is presented and discussed together with an 
interpretation involving simple topological concepts. Section 0 is devoted to some remarks 
and comments. 


II. MODELS AND NUMERICAL RESULTS 


To study the relationship between microscopic dynamics and equilibrium phase transi¬ 
tions we consider Hamiltonian systems of the standard type 


H[<p,n\ = ^ + , (!) 

i 

where (pi and Ti{ are canonically conjugated coordinates and momenta, i labels the sites of a 
d-dimensional cubic lattice, and V is an interaction potential. 

More precisely, we consider models which can be derived from the paradigm Hamiltonian 

H[(p] = jd d x j^vr 2 (x) + J^[V^(x)] 2 - i^ 2 (x) + ^ 4 (x)| • (2) 

where 7 r(x) = 6L[ip,ip]/Sip(^.) = <p(x) is the canonically conjugated momentum density of 
<£>(x), by discretizing it on a lattice. By means of the following substitutions 


3 Mx) _ + °e,) - gW _ 

a 

Jd d x —> a d , 


(3) 


we obtain 


= a d Y^ 


1 2 d 


i+eu 




fi=i 


1 99,^4 

2 m ^ + 4 ^ 


(4) 


where a is the lattice spacing, is the unit vector in the //-th direction of the lattice and 
ipi = (p(?t i). This system shows (at equilibrium) a continuous phase transition with nonzero 
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critical temperature corresponding to a spontaneous breaking of the discrete 0(1) — or Z 2 
— symmetry. 

We have also considered the vector versions of this lattice ip 4 model described by the 
Hamiltonian 


H [vM = ° d EE 


2™ 


+ ^I>f 

1 


i+e„ 




2 ™ 2 (<p?) 2 


A ^ 

+ zE 


I>f ) 2 


(5) 


where the index a runs from 1 to n. We have considered, in addition to n — 1, n — 2, which 
is the simplest vector case, and n = 4, which is the largest value of n that allowed for a 
complete numerical study with our computing resources. For n > 1 the broken symmetry 
is a continuous one (the potentials are respectively invariant under planar rotations, 0(2), 
and under the action of the 0(4) group). Because of the Mermin-Wagner theorem, the 
interactions being of short range, the 0(2) and 0(4) models can have a second order phase 
transition only on three dimensional lattices. 

The Hamiltonian dynamics — and thus the related dynamical, thermodynamical and 
geometrical quantities — is studied by molecular dynamics simulations performed at several 
values of the energy density £ = E/N , which is the relevant physical parameter as long as 
our systems are in a microcanonical ensemble 1 . 


A. Numerical study of dynamics and thermodynamics 


The canonical equations of motion 




d H 


■OL 


d H 
d'P? 


( 6 ) 

(7) 


1 The qualitative features of the results are not affected if we consider the temperature (average 
kinetic energy per degree of freedom) as the physical parameter. 
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yield 


with 


Vi = O 


■a 

7T; 


= ^ 2 > 

m=i 


a 

i+e M 


+ Vf-eJ + ^f-C7K 


IV? 


( 8 ) 


A = Ja d ~ 2 

B = m 2 a d - 2Ja d ~ 2 d (9) 

C = Xa d , 


and ||ipf || 2 = 'Yh a {Vi) 2 - I 11 order to guarantee a faithful numerical representation of a 
Hamiltonian flow, it is necessary that the algorithm updates the canonical coordinates, 
[<pf(nAt), (n^nAt)] —■» [<p“((n + l)At), (vr"((n + 1)At)], by means of a canonical, i.e. sym- 
plectic, transform. Symplectic algorithms ensure the conservation of Poincare geometric 
invariants, and, in particular, of phase space volumes and energy conservation. We used 
a very efficient and precise third order symplectic algorithm recently proposed 0 , keep¬ 
ing the fluctuations of relative energy at A E/E ~ 1CP 9 . All the simulations have been 
performed using words of 64 bits. We have always chosen random initial conditions at 
equipartition among momenta in order to consider phase space trajectories stemming from 
initial conditions that belong to the support of an equilibrium measure. 

Along the phase space trajectories — worked out numerically — the time averages of 
any observable A is computed as 


cItA[it(t), < p ( r )] . ( 10 ) 

By means of such averages both dynamical and thermodynamical properties of the system 
under investigation can be determined. 

One of the most relevant properties of the dynamics is its degree of instability, because 
it is related to the efficiency of phase mixing. Let us remember that the the strength of 
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dynamical instability, i.e. of chaos , is measured by the largest Lyapounov exponent Ai. If 
we denote by A4 the phase space of the system and by X a vector field on it, such that 

x 1 — X i (x 1 ... x N ) (11) 


are the equations of motion, a complete integral of this dynamical system defines a one- 
parameter group of diffeomorphisms of AA, that is (j) 1 : Ai —» AA. Denote by 

? = TMm k ( 12 ) 


the tangent dynamics equation, i.e. the realization of the mapping dcjfi : T X AA —> T^t^AA, 
where [J /| is the Jacobian matrix of [A"*], then the largest Lyapounov exponent Ai is defined 
by 


A, = .im bn »«*>» 


w001 lie(0)|| 


(13) 


and, by setting A[x(t),£(t)] = £ T J[x(t)] £/ = £ £/£ £ = ln(£ T £), this can be for¬ 

mally expressed as a time average 


Al = , lim V+ / dTA i. x ( r )it(. r )} ■ 

t xx> zt /n 


In practice, as we deal with standard Hamiltonians, the tangent dynamics 
written in the form 


£15 

dt 2 


+ 


d 2 V 


dWidyi ) " q 


e = o 


(14) 
can be 

(15) 


which, integrated along any numerical trajectory of Eqs.(|8]), makes possible the estimate of 
Ai from 

, a r 


^1 (k/v) — 


AfAt 


J2 ln 


n =1 


U(t n )\\ 

IK(*n-l)ll 


(16) 


where {£*} = ({£)}, {£*}), fp(i) = [C q {t + At) - fj(t - At)]/2At, and t n = nAt (At is some 
time interval). The average is extended up to a final time tj^ such that Ai(tjv') has attained 
a bona fide asymptotic value. 
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For what concerns thermodynamic observables, temperature — the basic quantity — is 
determined through the time average of kinetic energy per degree of freedom 



where N is the number of lattice sites and t is the total time during which a phase space 
trajectory is followed. This quantity shows a fast convergence in time and is expected to 
differ from its canonical counterpart by a O(j^) correction. 

Besides the bifurcation of the order parameter — (cp) — at some critical value of the 
temperature, a second order phase transition is signaled by a singular temperature depen¬ 
dence of the specific heat, and therefore the microcanonical computation of the constant 
volume specific heat Cy deserves special care. An efficient numerical method to compute 
Cy is devised by inverting a general formula relating canonical and microcanonical averages 
of the squared fluctuations of a generic observable [|19j by applying it to the fluctuations of 
kinetic energy 


SK2 =™ 2 -f v (w 


(18) 


where Cy = ( dE/dT ); overbar and hat stand for microcanonical and canonical averages 
respectively. 


The quantity 5K 2 can be easily computed along the numerical trajectories, whereas the 
analytic expressions K = K = N/2/3, 5K 2 = N/(2f3 2 ) are readily found. By inverting 
the equation above one immediately finds a formula for a microcanonical estimate of the 
canonical specific heat 


Cy = 


Nn 1 

2 1 - (Nn/2)[(6K 2 )/(K 2 )} 


(19) 


which requires the numerical computation of time averages of kinetic energy and of its 
squared fluctuations; Nn is the total number of degrees of freedom. 


B. Dynamical evidence of the phase transition 
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1. Detecting the transition: Binder cumulants 


In the canonical ensemble, a phase transition may show np only in the thermodynamic 
limit. As long as N is finite, all the thermodynamic quantities are regular functions of the 
temperature, and ergodicity and symmetry are not broken. Nevertheless, some marks of 
the transition show up neatly also in a finite system. The specific heat does not diverge, 
but exhibits a peak — whose height grows with the size of the system — at a temperature 
T^' V (N). In principle the order parameter is expected to vanish on the whole temperature 
range for any finite value of N, though in practice, e.g., in a canonical MonteCarlo simulation 
where the length of the sampling of yj is necessarily finite, the system is trapped in one of 
the two phases for a “time” which grows exponentially with N |J, and thus a fictitious 
symmetry breaking is observed at a temperature T£{N). This temperature, in general, does 
not coincide with T^ V (N), even if 


lim T? V (N) = lim T*(N) = T c °° . 


( 20 ) 


In the microcanonical ensemble ergodicity breaking may occur also at finite N, hence we 
can expect that a “true” critical energy exists also at finite N. No rigorous theoretical result 
is at our disposal regarding this aspect. Nevertheless, on the basis of asymptotic equiva¬ 
lence of statistical ensembles, the behaviour of microcanonical thermodynamic functions is 
reasonably expected to be similar to the canonical case, at least as N is sufficiently large. 
Indeed this is what is observed, as we shall see in the following. In particular, we expect the 
specific heat to exhibit a peak at a critical energy density which is a function of N. 

In the framework of the statistical theory of critical phenomena, by means of the finite- 
size scaling analysis |22| , |2Tf the critical properties of the infinite system are inferred from the 
values of the thermodynamic observables in finite samples of different sizes. In particular it 
is possible to locate the critical point by means of the so-called Binder cumulants 0. The 
Binder cumulant g that we have computed for our systems is defined as 

(</} 


9 = 1 


3(^ 2 ) 2 ’ 


( 21 ) 
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where 


^ 2 "> = ((!»« 

<»>)<. = ' 

i 

In the disordered phase the probability distribution of the order parameter will be nearly 
Gaussian with zero mean, hence g ~ 0. At variance to this, at zero temperature (or energy), 
when ipi = (po with no fluctuations, g = 2/3. At different sizes of the system, g will decay 
following different patterns <7 (IV, T) from 2/3 to 0 at increasing temperature. The remarkable 
fact is that the value of g at T c °° is independent of N, provided N is large enough for the 
scaling regime to set in; hence the critical point can be located by simply looking at the 
intersection of the different curves g(N, T) for different values of N. In principle, two different 
sizes are sufficient to locate the transitions; in practice, owing to the unavoidable numerical 
errors which affect g, it is necessary to consider at least three values of N. Moreover, the 
value of g at the critical point, usually referred to as g*, is a universal quantity, like the critical 
exponents; for a simple proof see e.g. Ref. [2T]. The importance of the Binder cumulant 
method is not only that it allows to easily locate the critical temperature, without the need of 
an extrapolation of the asymptotic behaviour of the fictitious finite-A" critical temperatures, 
but also that such an estimate of T c °° is independent of the other thermodynamic observables 
like (ip) or Cy, and this is obviously a great advantage in determining the actual critical 
behaviour — in particular the critical exponents. Moreover, one can regard the existence of 
a crossing of different curves g(N, T) as a “proof” for the existence of a phase transition in 
the system under investigation. This may be useful in various cases where the presence of 
singularities in the thermodynamic functions or the existence of a nonzero order parameter 
are difficult to observe (e.g. this is the case of spin glasses [[23|j). 

The theory behind the Binder cumulant method is totally internal to canonical statistical 
mechanics: to our knowledge no extension of this theory to the microcanonical ensemble 
exists. Nevertheless we will adopt the pragmatic point of view of assuming its validity as a 
numerical tool also in our dynamical simulations, and our operative definition of the critical 
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energy density e£° will be the intersection point of the curves g(N,e ) at different N. The 
consistency of the method will be checked a posteriori. In the following, unless explicitly 
stated otherwise, e c and T c will denote respectively and T c °°. 

The results for g(N, e) at different sizes for the if 1 lattice models are shown in Figs.[L]a, b, c. 
The crossing of the various curves at s c ~ 31 for the 0(1) model is quite evident, and similarly 
at e c ~ 44 for the 0(2) model, and at e c ~ 56 for the 0(4) model. 

Such estimates of the critical energy densities are obviously far from being extremely 
accurate. However, we are mainly interested in showing that the dynamical phenomenology 
is actually consistent with the existence of an equilibrium phase transition at finite energy 
density, and the values of e c are needed to understand whether or not the singular (or, more 
generally, peculiar) behaviours of the observables — either thermodynamical, or strictly 
dynamical, or geometric ones — that we are going to study can be associated with the 
phase transition. 


2. Temperature 

The temperature of the p 4 systems, numerically determined according to Eq. (0 , is 
plotted in Fig|| as a function of the energy density e. Note that for all the models a change 
of the function T{e) is clearly evident at e = e c . 

By plotting the Binder cumulants vs. the temperature T, the critical values T c are 
obtained for all the models and are found in complete agreement with the outcomes of the 
T{e) curves. These values are: T c ~ 35 for the 0(1) model, T c ~ 25 for the 0(2) model, 
and T c ~ 16 for the 0(4) model. 


3. Specific heat 

The specific heat c v = Cy/Nn per degree of freedom of the p 4 models here considered, 
computed according to Eq. (0). is plotted vs. the temperature in Fig.^j. The asymptotic 
values of the specific heat in the limits T —» 0 and T —> oo are exactly known. In fact at 
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low energies the anharmonic terms in the Hamiltonian can be neglected, thus the system 
behaves as a collection af harmonic oscillators and c v —> 1 as T —> 0. In the high-energy 
limit the quadratic terms in the potential are negligible with respect to the quartic ones, 
whence c v —> l/2 + l/4 = 3/4asT —> oo. At intermediate energy densities, neat peaks 
show up whose positions are close to T c for each model respectively. The heights of the 
peaks are found to grow with N and to decrease with n. 

4■ Dynamical properties 

We have shown that the outcomes of the dynamical numerical simulations of the scalar 
and vector versions of the lattice (p A model are perfectly consistent with the expectations of 
the effects of a second order phase transition on a finite sample. As already motivated above, 
this first result is non-trivial. Up to now its content is that — for all practical purposes — a 
dynamical simulation is actually equivalent to a microcanonical one, so that, at sufficiently 
large N, the results are in natural agreement with canonical statistical mechanics. All 
these results concern time averages: the time variable, even if not eliminated from the very 
beginning as in the statistical approach, has been nonetheless integrated out in the averaging 
procedure. But we can also wonder what are the properties — if any — that are peculiar to 
the dynamics and that can be considered relevant to the description of the phase transition 
itself. Moreover we have already noticed that the phenomenon of ergodicity breaking has a 
deep dynamical origin, therefore we can try to understand what features are associated to 
a Hamiltonian ergodicity breaking. 

The lattice <p 4 models under investigation are nonintegrable dynamical systems. In the 
two limits £ —> 0 and £ —> oo, these systems become integrable. The two integrable limits 
are respectively those of a system of coupled harmonic oscillators and of a system of inde¬ 
pendent quartic oscillators. The dynamics is always chaotic over the whole energy range. 
Nevertheless, in analogy to other nonlinear oscillator systems, by varying the energy we 
expect that qualitatively different dynamical regimes will be found, characterized by a tran- 
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sition between different behaviours of the largest Lyapounov exponent Ai as a function of 
energy density or, equivalently, temperature. This phenomenon is attributed to a dynamical 
transition between weak and strong chaos, it is known as Strong Stochasticity Threshold 
(SST) and is discussed in refs. p4| , |25| . In particular, the following questions naturally arise. 
Is there any peculiar behaviour of the Lyapounov exponent in correspondence with the 
phase transition? Is there a transition between strongly and weakly chaotic regimes also in 
these models, and, in the affirmative case, is there any relationship between these different 
dynamical regimes and the thermodynamic phases? 

We must say from the very beginning that there are not yet conclusive answers to these 
questions. The study of a possible relation between chaos and phase transitions is a very 
recent issue |26|, and the results so far obtained and reported in the literature range from the 
claim of the discovery of a “universal” divergence in Ai near criticality in a class of models 
describing clusters of particles 0 , to the observation that the Lyapounov exponent attains 
its minimum in correspondence with the phase transition in Ising-like coupled map lattices 
and to the apparent insensitivity to the liquid-solid phase transition of the Lyapounov 
spectra of hard-sphere and Lennard- Jones systems [ jZ£f |. 

Our simulation results are plotted in Figs.[| and |5|. The 0(1) case has been studied more 
extensively than the others because of practical reasons of computational effort (for example 
single runs for the 0(4) model usually required at least two weeks of CPU time on a fast 
HP 9000/735 computer). 

The first numerical evidence is that in presence of a second order phase transition a rather 
sharp and “cuspy” transition between different behaviours of Ai(T) is found at T c (where 
the critical values T c are those determined by means of Binder cumulants). Moreover, the 
qualitative behaviour of Ai(T) appears very different in the thermodynamically ordered and 
disordered regions respectively. In fact, in the former region Ai rapidly increases with T, 
whereas in the latter region Ai(T) displays an almost flat pattern above T c (note that Ai(T) is 
expected to change again at very large T because the dynamics is asymptotically integrable 
in the limit T —> oo; this effect has been numerically checked at very high temperatures and 
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is clearly evident in Fig.|4| — 0(1) case — at T/T c ~ 10 4 ). This suggests that the phase 
transition has a dynamical counterpart in a passage from a weakly to a strongly chaotic 
regime. 

It is remarkable that the shape of Xi(T) is significantly different in presence or in absence 
of a second order phase transition. In fact, in the case of one dimensional lattices with short 
range interactions — where no phase transition is present — Ai(T) has a very smooth pattern 
(see Ref. ||). This fact has been checked more specifically for the <^ 4 model by computing 
Ai(T) for the 0(2) symmetry case on a two-dimensional lattice; as a consequence of the 
Merrnin-Wagner theorem, here a second order phase transition is forbidden and in fact this 
model undergoes an infinite order phase transition (Kosterlitz-Thouless-Berezinsky). The 
shape of Xi(T) again displays a major change so that the low and high temperature regimes 
are very different. However the transition between these two regimes is now smooth [|30| . 


It is worth to emphasize that the average of a local property of microscopic dynamics 
- the average instability measured by Ai — is sensitive to a collective phenomenon like a 
second order phase transition. 

It could be argued that in the critical region almost any “honest” observable will show a 
peculiar behaviour and that this reflects the tendency of the statistical measure to become 
singular at the transition point, regardless of the ensemble chosen. In the framework of 
equilibrium statistical mechanics this is certainly true, because Gibbs measure is the funda¬ 
mental mathematical object upon which everything relics. In the thermodynamic limit also 
the microcanonical measure, which is the invariant measure of the microscopic Hamiltonian 
dynamics, will have to become singular. However the microcanonical measure is not the 
ultimate mathematical entity that can be considered, so that the Hamiltonian dynamics 
approach gives meaning to the question of the possible existence of a more fundamental 
phenomenon at the very ground of a phase transition. 

Lyapounov exponents provide the necessary link to such unexplored land. The details 
on this point are given in the next Section, where we recall how the geometrization of 
Hamiltonian dynamics proceeds in the language of Riemannian geometry, and how average 
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geometric properties of some suitable manifold directly influence the average dynamical 
instability quantified by Ai. 


III. GEOMETRY OF DYNAMICS AND THE PHASE TRANSITION 

Let us here sketch the main points of the Riemannian theory of chaos in physical systems, 
details can be found in ref.s |3|-|8|. 


A. Riemannian geometrization of newtonian dynamics 


The trajectories of a dynamical system described by the Lagrangian function 

£(q> q) = ^(q )q l q k - V (q) (22) 

are geodesics of the configuration space endowed with a proper Riemannian manifold struc¬ 
ture described by the metric tensor 


Satq) = 2 [E - V(q)]o»(q) . 


(23) 


This metric is known as Jacobi metric and is defined in the region of the configuration space 
where E > R(q). In local coordinates, the geodesic equations on a Riemannian manifold 
are given by 


d ¥ ri dq^dq^_ 
ds 2 ]k ds ds 


(24) 


where s is the proper time and T l - k are the Christoffel coefficients of the Levi-Civita 
connection associated with , i.e. T* fc = ■^8 im (djW8km + dkWS m j — d m W5jk), where 
W = E — V( q); proper time and physical time are related by ds 2 = 2 W 2 dt 2 . By direct 
computation, using g & — (E — V (q)) it can be easily verified that the geodesic equations 
yield 


d 2 q i _ _dV 
dt 2 dq i 


( 25 ) 
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i.e. Newton’s equations associated with the Lagrangian (^). These equations can be also 
derived as geodesics of a manifold consisting of an enlarged configuration spacetime Mxl 2 , 
with local coordinates (q°, q 1 ,..., q\ .. ., q N , q N+1 ). To such purpose this space is endowed 
with a non-degenerate pseudo-Riemannian metric, first introduced by Eisenhart [[3TJ, whose 
arc-length is 

ds 2 = g, JiV dq^dq v = a tJ dq l dqi — 2V (q )(dq 0 ) 2 + 2 dq°dq N+1 , (26) 

called Eisenhart metric. The natural motions are obtained as the canonical projection of 
the geodesics of (M x M 2 ,^) on the configuration space-time, i : M x R 2 h M x R. 
Within the totality of geodesics only those whose arclcngth is positive-definite and is given 
by ds 2 = cfdt 2 correspond to natural motions, what is equivalent to requiring the condition 
q N+1 = \c\t + c\ — f* Ldr for the extra-coordinate q N+1 PH; c\ and C 2 are real arbitrary 
constants. 


B. Curvature and instability of geodesic motions 


There is an important relation between the curvature of a manifold and the stability of 
its geodesics. It is described by the Jacobi - Levi-Civita (JLC) equation for the geodesic 
separation vector field J(s). 

The evolution of J contains the whole information on the stability — or instability — of 
a given reference geodesic y(s): in fact if \J\ grows exponentially then the geodesic will be 
unstable in the Lyapounov sense, otherwise it will be stable. It is remarkable that such an 
evolution is completely determined by the Riemann curvature tensor R l - kl according to the 
JLC equation 


ds 2 


o i d d j jk d Q l 

K > kl ds J ds 


= 0 


where y- is the covariant derivative. 

as 


(27) 


In the large N case, under suitable hypotheses []7|,|8|], it is possible to derive a scalar 
effective stability equation. Briefly, among the others, the main assumptions are: i ) that 
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the ambient manifold is almost-isotropic, this essentially means that - after some suitable 
coarse-graining - the ambient manifold would look like a constant curvature manifold; ii) 
that the curvature felt along an unstable geodesic can be reasonably modeled by a gaussian 
stochastic process. The final result is || 

^ + {S 2 k R )l /2 ri(s)ip = 0 , (28) 

where ip denotes any of the components of J in Eq. (|27|) because now all of them obey 
the same effective equation of motion; (/cr)^ = where Kr is the Ricci curvature 

of the ambient manifold: Kr = Rikq l q k and R,j. = i?b fc ; (-) M stands for microcanonical 

average, and {S 2 kR) tl is a shorthand for jA^(S 2 Kr) 11 , the mean square fluctuation of the Ricci 
curvature; 77 (s) is a gaussian 5-correlated random process of zero mean and unit variance. 

Equation (|28|) is a scalar equation which, independently of the knowledge of dynamics , 
provides a measure of the average degree of instability of the dynamics through the growth- 
rate of ip(s). The peculiar properties of a given Hamiltonian system enter Eq. (^8|) through 
the global geometric properties (/cr)^ and (<5 2 /cr of the ambient Riemannian manifold. 
Moreover (/cr)^ and {5 2 kR) il are functions of the energy E of the system — and of the 
energy density e = E/N as well, which is the relevant quantity at N —» 00 — so that from 
(^H) we can obtain the energy dependence of the geometric instability exponent. 

Equation (^) is of the form 

44 + ^)^ = 0 (29) 

ds z 

representing a stochastic oscillator where the squared frequency fi(s) is a stochastic process; 
the derivation of this equation does not depend on a particular choice of the metric. For 
Hamiltonian systems with a diagonal kinetic energy matrix, i.e. a tJ = S tJ , by choosing as 
ambient manifold for the geometrization of dynamics the enlarged configuration space-time 
equipped with Eisenhart metric (^6|), it is found that the only non-vanishing component of 
the Ricci tensor is i?oo = AH, thus Ricci curvature is a function of the coordinates q 1 only 
and one has kR^q) = AV/N. Using dt 2 = ds 2 , the stochastic oscillator equation fl2§| ) can be 
written 
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^ + n(t)* = o, 


( 30 ) 


where mean and variance of Q(t) are given by 

% = {k R )„ = i(AV}„ , (31) 

<4 = <^r)^ = ^ «(AV') 2 )„ - (A V)l) . (32) 


The process Q(t) is specihed by flo, and its time correlation function Tn^i,^)- We 
consider a stationary and ^-correlated process Q(t) with r f1 (f 1 ,t 2 ) = roy) <5(|£ 2 — £i|) , where 
r is a characteristic time scale of the process. At present the evaluation of this time scale 
is still a rather delicate point, where some arbitrariness enters the theory. In Ref. || these 
two time scales are defined by 


T, = /*\ 


7r 


t ds / 2VV 0 + an 
that are combined to give r as follows 


t 2 = { -T 


dt \ nl /2 


ds ctq 


(33) 


r 


-l 


2 (r, 1 + t 2 ') 


(34) 


As we shall see below, at low temperatures this formula seems to predict a satisfactory 
temperature-dependence of r, in fact, by adjusting a constant factor that multiplies Ti, the 
theoretical prediction of Ai(T) is in very good agreement with numerical computations. At 
high temperatures we have to take care of the fact that flo and an are both increasing 
functions of T, even though the system approaches an integrable limit. 

Whenever Q(t) in Eq. ([3(]) has a non-vanishing stochastic component, the solution -0(f) 
is exponentially growing on the average [Q. Our estimate for the (largest) Lyapounov 
exponent is then given by the growth-rate of ||(^,'0)(f)|| 2 according to the definition 


A, = l im Ilog.^M+A*) 


(35) 


t^oo2t t /; 2 ( 0 ) + V > 2 ( 0 ) 

The ratio ('0 2 (f)+'0 2 (f))/(?/> 2 (O) + '0 2 (O)) is computed by means of a technique developed 
by Van Kampen, summarized in Ref. R , that yields 
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Ai(f2 0 , ctq,t) 


= 2 IA 


4^o 

~3A 


1/3 


A = 


2 <r a T + 


4 a 


+ t 2 ^) 2 


( 36 ) 


The quantities h2 0 , cr^ and r can be computed as static, i.e., microcanonical averages. There¬ 
fore Eq. ( |36| ) gives an analytic, though approximate, formula for the largest Lyapounov 
exponent independently of the numerical integration of the dynamics and of the tangent 
dynamics. 


C. Geometric signatures of the phase transition 

As already noted above, on the one hand the largest Lyapounov exponent is sensitive 
to the phase transition, on the other hand, in the Ricmannian description of chaos, Ai is 
intimately related to the average curvature properties of the “mechanical manifolds”. These 
quantities are computed as integrals on manifolds just like other statistical quantities of 
thermodynamic kind. This means that by means of statistical-mechanical-like computations 
we can obtain non-trivial informations about dynamics. Hence the following questions: is 
there any peculiarity in the geometric properties associated with the dynamics of systems 
which — as statistical systems in thermal equilibrium — exhibit a phase transition? And in 
particular, do the curvature fluctuations show any noticeable behaviour in correspondence 
with the phase transition itself? 

1. Results of the computations 

Let us now report on the results of the the computation of the geometric properties of 
the “mechanical” manifolds sampled by the numerical geodesics. 

For the </d models, the Ricci curvature per degree of freedom along a geodesic of (M x 
M 2 , < 7 e) is given by 
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k R = a7" E E 55TT2 = 2J(i - “ 2 + A Tvr^EE(»’f) 2 • 


(37) 


a=l i v i 7 a=l i 

High and low temperature behaviours of this quantity can be easily derived. In the limit 
T —> 0 we can replace, at any site i of the lattice, ||^"|| 2 = with the constant 

value tpl = m 2 / A; this value is obtained by minimizing the potential part of the Hamiltonian 
(H). Hence, for a generic 0(n ) case, we have 


lim kr> = —(Jdn + m 2 ) 
T—>o n 


(38) 


and with the values we chose for the constants — J — 1, m 2 = 2, d = 3, A = 0.1 - 
it is kn = 10 in the 0(1) case, kji — 8 in the 0(2) case, and kn = 7 in the 0(4) case 
respectively. These values perfectly check with our numerical findings as is shown by Fig.|6] 
where k(T) = (( kR)t(T ) — 2Jd)/((kn) t (T = 0) — 2 Jd) is synoptically displayed for all the 
models; the average (•)* is defined in Eq. (|10|). At low temperature (/c#)t(T) only slightly 
deviates from its limiting zero-temperature value, as is shown by Fig.[7j; this fact is intuitively 
interpreted as a sign of a weakly chaotic dynamics. 

Also in the opposite limit, T —> oo, these systems are again integrable. In fact at 
increasing temperature the variables 99 " become larger and larger so that the Hamiltonian (^|) 
describes a collection of quartic oscillators that are less and less perturbed by the quadratic 
coupling term. In this limit the canonical partition function is factored in terms of functions 
of the following form 

[ dxx u e~ axa = -r ( a - ( " +1)/ “ , (39) 

JO OL \ OL J 

with v — 0, and where F(x) is the Euler gamma. Hence the canonical average of any even 
power v of the field is 

•/sav* r(^) 


MY) = 


(40) 


4/ r (j) 

and vanishes for any odd power v. 

^Frorn Eqs.(l37|) and (j4(J) we find a canonical estimate of ( k ^which differs from the 
microcanonical one by 0(A) terms 
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(41) 




2 Jd 


m 2 + 


2(n + 2)r(^)VX 

r (l) 


VT +0(-) 


This prediction is compared to the numerically computed values of (k R ) t {T) in Fig.|7|: at 
very high temperature the agreement is very good. 

In order to compute the average curvature fluctuations, (S 2 k R ) ll , we first notice that 

(f-kn) = (kl) - (k R ) 2 = A 2 INI 2 ) ''j - (£<|M| 2 >) I (42) 

and, as in the large T limit we consider all the decoupled, we find 


(i 2 fc a )=A 2 (n + 2) 2 {((^, ) r)-((^, ) ) 2 ) 2 } 


(43) 


where <p>%. denotes any representative of the now independent degrees of freedom. The 
Gibbsian, canonical, average in the T —> oo limit is now easily found to be 


(■ S 2 k R ) 


G 


( £(5/4) fr(3/4) 


(n + 2) 2 4A T . 


(44) 


\T(l/4) [T(1/4)J 

In order to compare the predictions of Eq. (f|4|) with our numerical results, and also in order 
to use it in the analytic prediction of the Lyapounov exponent, we have to take into account 


the correction that relates canonical and microcanonical averages [p2|] that now reads 


\ l&h \G P 2 ( d ( k R)X 

(3 k R )„ = (6 k R ) 


(45) 


The high temperature partition function Z is obtained by raising to the iVri-th power 
the integral f d(/?exp[— f3{\/A)(p A ] ~ Then, using F = — (l/Nn/3) InZ and Cy = 

—T(d 2 F/dT 2 ), we find c v —> 1/4. By the way, this is in very good agreement with our 
numerical results for the high temperature values of c v ; this is somehow less clear in the 
0(4) case because c v was computed only in the transition region. From Eqs. (|^) and 
we can now obtain the final result 


(S 2 k R ), 


f £(5/4) 

\r(i/4) 


r(3/4) 

r(i/4) 


(n + 2) 2 4A T . 


(46) 
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In Fig.|8] we report the temperature dependence of the time average of the Ricci curvature 
fluctuations, ctq{T) = (5 2 k R ) t . In Fig.[|we also give a comparison of ctq with the prediction 
of Eq.(f46|) for the 0(1) model. 

The common feature of the three models is that a cusp-like (singular) behaviour of the 
curvature fluctuations is observed in correspondence with the phase transition. 

Moreover, curvature fluctuations display very smooth energy density dependence, or 
temperature dependence as well, in those systems where no finite order phase transition is 
present (see ref. ||). In Fig.[HJwe report also an(T) in the case of 2-d <p 4 model with 0(2) 
symmetry; in this case a second order phase transition is forbidden and actually the system 
undergoes a Kosterlitz-Thouless phase transition. The cusp-like behaviour of curvature 
fluctuations has now disappeared and <7 q(T) is a monotonically increasing function of T; 
visibly, something still happens at the transition point (T c ~ 1.5) so that this case appears 
to be “intermediate” between no phase transition at all and a second order phase transition. 


Similar results have been found for planar 2-d and 3-d classical Heisenberg models |30|J33 
and in a preliminary investigation of the dual (gauge) version of the Ising model in 3-d pi 


the cusp-like behaviour of the curvature fluctuations always shows up when a second order 
phase transition is present, and the singular point is located at the critical temperature, 
within the numerical accuracy. 

In the light of the Riemannian description of Hamiltonian chaos given above, we under¬ 
stand why the temperature dependence of the largest Lyapounov exponent Ai is so peculiar 
near and at the critical temperature (see Figs.[| and[5|): Ai(T) reflects the “cuspy” pattern 
of oq(T) near T c . In the next Subsection we make a conjecture about the deep meaning of 
these singular behaviours shown by Ai(T) and ctq{T). 

As the invariant measure for an autonomous Hamiltonian flow is the micro canonical mea¬ 
sure on the constant energy surfaces of phase space, our numerical computations of (k R ) t 
and of (5 2 k R ) t are good estimates of the quantities Ho(F) and an(T), i.e. microcanonical 
averages, that enter Eqs. i fTS| ) and (^). The analytic computation of Ai(T) by means of 
these formulae yields an unsatisfactory result that overestimates Ai(T) at low temperatures 
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(though the temperature dependence is correct) and that steeply increases at high tem¬ 
peratures instead of saturating (before decreasing again at extremely high T). The high 
temperature result appears particularly bad, however this is only due to the asymptotic 
growth with T of both (kp) and ( S 2 kp ) — given by Eqs. ( [41]) and ( |46| ) and confirmed 
numerically — which has no special meaning for dynamical instability. The estimate of 
the decorrelation time-scale of curvature fluctuations along the geodesics is still somehow 
rudimental in the above outlined Riemannian framework, therefore one expects that some 
improvement is needed on this point. As a matter of fact, it is possible to substantially 
improve the theoretical predictions by simply multiplying the decorrelation time scale r of 
Eq.(|34|) by a constant factor which is model dependent and different below and above T c . 
Moreover, at high temperatures, in computing T\ and r 2 given in Eq. (|33|) , we have subtracted 
to (T) and ctq(T) their respective asymptotic behaviours given by Eqs. (|4]) and (fl6|). 
The analytic predictions for Ai (T) are now in very good agreement with numeric results with 
the exception of the critical region, where something is apparently still lacking. The results 
are reported in Figs. Q.0 and 0 where it is well evident that the best agreement between 
theory and numerical experiments is obtained in the 0(1) case; a very good agreement is 
still present at low temperatures for the 0(2) model and it becomes poorer in the 0(4) 
model. The comparison at T > T c suffers - in the cases of 0(2) and 0(4) - of a restricted 
range of temperature values (we focused our attention only to the transition region because 
of the already mentioned problems) where subtracting to f! 0 (T) and ctq(T) their asymptotic 
values is less meaningful. 

However it is not out of place to remind that the theoretical computation of Lyapounov 
exponents is not a routine task at all, and that the approach reported here is at present the 
only theoretical method available to cope with the computation of Ai. What is important 
here is that with some simple and reasonable adjustment the above sketched analysis still 
applies and yields good results. Refinements of the geometrical theory of chaos are beyond 
the aim of the present work, rather we are interested in using it as it is at present to get a 
hold of the deep origin of the peculiarities of the dynamics at a phase transition. 
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D. A topological conjecture 


We shall now try to grasp the possible significance of the above reported cusp-like, thus 
possibly singular, behaviour of the curvature fluctuations at the transition point for the <p 4 
lattice systems. As a first step toward this goal, we shall try to reproduce such a peculiar 
behaviour of curvature fluctuations in abstract geometric models. A preliminary step in this 
direction was already presented in Ref. |3(J, applied to the case of planar spin models. 

The choice of a geometric toy model stems from the following considerations. Weakly 
and strongly chaotic geodesic flows can “live” on homologically trivial manifolds, i.e., on 
manifolds that are diffeomorphic to an Wsphere, in other words, non-trivial topology is 
not necessary to make chaos, conversely, a sudden topological change in a family of mani¬ 
folds can abruptly affect their geometric properties and the degree of chaoticity of geodesic 
flows. Therefore, let us consider — for instance — the two families of surfaces of revolution 
immersed in R 3 defined as follows: 


T e = (f e (u) cosv, f e (u) sin v,u) , (47a) 

Q e = {u cos v, u sin v, f £ (uj) , (47b) 

where 

fe{u) = ±V£ + U 2 -u 4 , £ G [£ m i n , +Oo) , (48) 

and e min = — 4. Some members of the two families are depicted in Fig. [14]. In both cases 
there exists a critical value of the parameter e, e c = 0, corresponding to a change in the 
topology of the surfaces. In particular the manifolds T e are diffeomorphic to a torus T 2 for 
e < 0 and to a sphere § 2 for e > 0. In the other case, one has instead a change in the 
number of connected components: the manifolds Q e are diffeomorphic to two spheres for 
e < 0 and to one sphere for £ > 0. Computing the Euler-Poincare characteristic x one 
finds y(7T) = 0 if £ < 0, and x(^e) = 2 otherwise, while x(Ge) equals either 4 or 2 when 
£ is respectively negative or positive. Let us now compute the e-dependence of the average 
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curvature properties of these surfaces as e —> e c . Let M belong to one of the two families 
under investigation. The gaussian curvature K is given by |35| 

x'(x"y' — x'y") 


K = 


y[x' 2 + y' 2 ) 2 


(49) 


where the functions x(u) and y{u ) represent the coefficients of the general form M(u,v ) = 
(y(u) cos v,y(u) sin v,x(u)) of parametrized surfaces of revolution, and the prime denotes 
differentiation with respect to u. Now the fluctuations of K are computed as follows 

2 


a 2 = ( K 2 ) - (K) 2 = A 


-l 


K 2 dS- A 


IM 


1-1 


IM 


KdSj , 


(50) 


where A is the area of M and dS is the invariant surface element. Both families of surfaces, 
in spite of having very different curvature properties on the average^, exhibit a singular 
behaviour in the curvature fluctuation a as e —> e c , as shown in Fig. |1^. 

These results suggest — at a heuristic level — that, from the point of view of the geomet¬ 
ric description of the dynamics, a phase transition might correspond to a topology change in 
the manifold underlying the motion. The relevance of topological concepts for the theory 
of phase transitions has been already emphasized (see Ref. |H^|) though in a more abstract 
context. Here we suggest that topological aspects of phase transitions might also concern the 
manifolds that are “just behind” dynamics, and not only those deep mathematical objects 
that are involved in Ref. |56). In our opinion this subject deserves further investigation to go 
beyond the heuristic level. In fact the study of dynamics and of its geometric and topologic 
counterparts could eventually lead to a better understanding of the nature of ergodicity 
breaking and thus of the very nature of phase transitions. 


2 For instance, (K)(e) = 0 in the T e case as e < 0, while the same average curvature is positive 
and diverging as e —> 0 for Q e . 
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IV. CONCLUDING REMARKS 


Let us now summarize the main points of the present work and comment about their 
meaning. 

By studying some classical lattice ip 4 models that undergo second order or Kosterlitz- 
Thouless phase transitions, it has been found that the natural microscopic dynamics - 
derived from the Hamiltonian functions of these systems — clearly reveals the presence of 
the phase transition. The invariant measure of Hamiltonian dynamics is the microcanonical 
measure, equivalent — in the thermodynamic limit — to the canonical measure which is 
sampled by usual MonteCarlo algorithms. Therefore one could argue that it is not surprising 
that Hamiltonian dynamics yields the same results of a MonteCarlo stochastic dynamics. 
As a matter of fact using Hamiltonian dynamics just to sample the microcanonical measure 
would not be so interesting, whereas the important point raised by the present work is that 
Hamiltonian dynamics brings about new observables and a new framework to tackle phase 
transitions. Mainly Lyapounov exponents are the new observables intrinsic to the dynamics, 
and the differential-geometric treatment of dynamical instability is the new framework. 

As well as thermodynamic observables, dynamic and geometric observables are sensitive 
to a second order phase transition which can be recognized through their peculiar “non¬ 
smooth’' behaviours. The common wisdom on phase transitions suggests that non-smooth 
behaviours of any observable are expected near the transition point, as a consequence of the 
tendency of the measure to become singular at T = T c in the limit N —> oo. In the light 
of our results we suggest that a deeper explanation might be possible: a major topological 
change of the “mechanical” manifolds could be the common root of the peculiar behaviours 
of both dynamic and thermodynamic observables in presence of a phase transition. Here 
topology is meant in the sense of the De Rham’s cohomology. 

On a purely phenomenological ground it might be surprising that the largest Lyapounov 
exponent, which measures an average local property of the dynamics, is sensitive to a col¬ 
lective, and therefore global, phenomenon like a phase transition. In fluids, for example, 



it is well evident that molecular chaos has nothing to do with the macroscopic patterns of 
the velocity held. It is even possible to have chaotic motions of fluid droplets (Lagrangian 
chaos) in presence of regular Eulerian velocity fields (i.e. in laboratory reference frame). 

However, within the Riemannian framework outlined in the preceding Sections, Lya- 
pounov exponents appear tightly related to the geometry of the “mechanical” manifolds, 
and geometry dramatically changes in presence of a major change of topology. Thus our 
topological conjecture seems to naturally account for this - at first sight counterintuitive - 
sensitivity of the largest Lyapounov exponent to a macroscopic collective phenomenon. 

It is worth mentioning here that, at the best of our knowledge, there is only another 
framework where Lyapounov exponents can be, at least in principle, analytically computed. 
This is a field-theoretic framework, already mentioned in the Introduction |[T3|-[T5|, based on 
a path-integral formulation of classical mechanics, where Lyapounov exponents are seen as 
expectation values of suitable operators. There are many interesting points in this framework 
that could probably reveal a fertile relationship with the Riemannian geometric approach 
which is behind our present work. Let us mention some of them: ergodicity breaking, which 
- as we discussed in the Introduction - is a more general concept than symmetry breaking, 
in the field-theoretic context appears related to a supersymmetry breaking, moreover this 
supersymmetry breaking can occur also at finite N [|K|]; Lyapounov exponents turn out 
related to mathematical objects that have many analogies with definitions and concepts of 
the de Rham’s cohomology theory |jT5[ what might be useful in future investigations about 
the relation between Lyapounov exponents and topology at a phase transition. 

In conclusion we believe that the dynamic approach, besides the conceptual aspects 
mentioned above, could contribute to complement the standard approaches of statistical 
mechanics to the description of phase transitions, and it might hopefully be particularly 
helpful in those cases where these standard methods may encounter some difficulty, as is the 
case of disordered and frustrated systems, polymers in the continuum, lattice gauge theories 
where no symmetry-breaking transition occurs. 
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FIGURES 

FIG. 1. Binder cumulants g(N,e ) vs energy density e at different values N of the lattice sites. 
a) the 0(1) case; b) the 0(2) case; c) the 0(4) case. Open circles refer to JV = 4 x 4 x 4, full 
triangles refer to N = 6 x 6 x 6, and full circles refer to N = 8 x 8 x 8. 

FIG. 2. Temperature T (twice the average kinetic energy per particle) is plotted vs energy 
density s. Results of the 0(1), 0(2) and 0(4) models are represented by full circles, open circles 
and open triangles respectively. Temperatures and energy densities of each model are scaled by the 
corresponding critical values obtained by means of Binder cumulants. The dashed line is a guide 
to the eye. Lattice size: IV = 8x8x8 


FIG. 3. Specific heat per degree of freedom vs scaled temperature T/T c . cy = Cy/Nn and Cy 


is computed according to Eq.(19). Symbols: full circles for 0(1); open circles for 0(2) and open 
triangles for 0(4). Lattice size: IV = 8 x 8 x 8. 


FIG. 4. The largest Lyapounov exponent Ai is plotted vs temperature for the 0(1) model. A 
“non-smooth” feature at T = T c is well evident. Lattice size: IV = 8 x 8 x 8. 

FIG. 5. Synopsis of Ai(T) obtained for the 0(1) model (full circles), for the 0(2) model (open 
circles) and for the 0(4) model (open triangles). Lattice size: N = 8 x 8 x 8. 


FIG. 6. Reduced average Ricci curvature k = ((kn) — 2 Jd)/((kit(T = 0)) — 2 Jd) vs T/T c . Ricci 
curvature is so reduced in order to facilitate the comparison between the different models. Symbols: 
full circles for 0(1); open circles for 0(2) and open triangles for 0(4). Lattice size: N = 8x8x8. 


FIG. 7. 0(1) model. The average Ricci curvature (&#) is plotted vs T for a wide temperature 

range. The dashed horizontal line represents the integrable limit behaviour of (kn){T) predicted by 
Eq.(|3^) and actually attained at low temperature by the average Ricci curvature computed for the 


0(1) model according to Eq.(37). Solid line represents the high temperature asymptotic behaviour 
of ( kn)(T) predicted by Eq.(41), again pertaining an integrable limiting behaviour of the model. 


34 





FIG. 8. Average Ricci curvature fluctuations an vs T/T c . The “cuspy” behaviour is well evident 
at T ~ T c . From top to bottom: 0(4), 0(2) and 0(1) results. The cusp appear to soften at 
increasing dimension n of the symmetry group O(n). 

FIG. 9. Average Ricci curvature fluctuations an vs T for the 0(1) model reported for a wide 
range of temperature. Solid line represents the high temperature asymptotic value given by Eq.(f46|). 

FIG. 10. Average Ricci curvature fluctuations an vs temperature for the 0(2) model on a square 
lattice (d = 2) of N = 30 x 30 sites. The cusp is now absent and an CO is a monotonously increasing 
function. Around T ~ 1.5, on the basis of the temperature behaviour of other observables, the 
system is supposed to undergo a Kosterlitz-Thouless phase transition and, correspondingly, we can 
observe a change in the shape of <jq(T). Here J = 1, A = 4 and m 2 = 10. 

FIG. 11. The numerical largest Lyapounov exponent Ai (open circles) is plotted vs T for the 
0(1) model and is compared to the analytic prediction of Eq.(|36|) (full circles). The vertical solid 
line marks the transition temperature. The correlation time scale r is given by Eq.([34|); r is rescaled 
by a constant factor equal to 0.65 at T < T c , and by a factor 1.1 at T > T c . 

FIG. 12. The numerical largest Lyapounov exponent Ai (open circles) is plotted vs T for the 
0(2) model and is compared to the analytic prediction of Eq.([36|) (full circles). The vertical dashed 
line marks the transition temperature. Here r is rescaled by a constant factor equal to 3 at T < T c , 
and by a factor 0.7 at T > T c . 

FIG. 13. The numerical largest Lyapounov exponent Ai (open circles) is plotted vs T for the 
0(4) model and is compared to the analytic prediction of Eq.(|36|) (full circles). The vertical solid 
line marks the transition temperature. Here r is rescaled by a constant factor equal to 5.5 at 
T < T c , and by a factor 0.6 at T > T c . 
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FIG. 14. Some representatives of the two families of surfaces J- t and Q e defined in Eqs. fl47a| ) and 
(|47b| ) respectively. Each family is divided into two subfamilies by the critical surface corresponding 
to e c = 0 (middle members in the picture). Members of the same subfamily are diffeomorphic, 
whereas the two subfamilies are not diffeomorphic between them. 

FIG. 15. The second moment of the gaussian curvature of the surfaces and Q e is plotted vs 
e. a is defined in Eq.(|50|), e is shifted by e m i n = 0.25 (see text) for graphical reasons, (a) refers to 
Q t and (b) refers to The cusps appear at e = 0 where the topological transition takes place. 
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